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• Abstract. We study the neutron-rich calcium isotopes performing shell-model calculations 

' with a realistic effective interaction. This is derived from the CD-Bonn nucleon-nucleon potential 

£h . renormalized by way of the V[ ow _k approach, considering 48 Ca as an inert core and including 

the neutron Ogg/2 orbital. We compare our results with experiment and with the results of a 
previous study where 40 Ca was assumed as a closed core and the standard lpO/ model space was 
, employed. The calculated spectroscopic properties are in both cases in very good agreement with 

' the available experimental data and enable a discussion on the role of the gg/2 single-particle 

i state in the heavy-mass Ca isotopes. 

d 

■ 

1. Introduction 

Heavy-mass calcium isotopes (N > 28) have been recently the subject of a renewed experimental 
and theoretical attention, since they lie far from the stability valley, thus allowing to explore the 
evolution of the shell structure when approaching the neutron drip line. 

This has resulted in a large number of experiments in this region [1-10] aiming to 
study the evolution of the single-particle (SP) orbitals. On the theoretical side, to describe 
the spectroscopic properties of these nuclei the shell model with various two-body effective 
interactions has been extensively used in recent years. The most frequently employed interactions 
are the KB3G [11], FPD6 [12], GXPF1 [13], and a new version of the latter dubbed GXPF1A 
[14]. The above interactions have all been derived by modifying the two-body matrix elements 
(TBME) of some realistic effective interactions in order to obtain a better description of the 
experimental data. 

During the last years, however, realistic shell-model effective interactions without any 
empirical modifications have proved to be able to describe with remarkable accuracy the 
spectroscopic properties of nuclei in various mass regions [15-18]. 

On these grounds, we have recently performed a shell-model study of the Ca isotopic chain 
[19] employing an effective interaction derived from the CD-Bonn potential renormalized by way 
of the Mow-k approach [20, 21]. In that work we have chosen a model space spanned by the four 
neutron SP levels of the fp shell outside the doubly-closed 40 Ca core. This is a standard choice 
in most of the studies on neutron-rich Ca isotopes, and has therefore allowed a direct comparison 
between our results and those of other shell-model calculations. 

In the present work, starting from the same Mow-k as m R- ei - [19], we have derived a shell- 
model effective interaction for a model space spanned by the neutron SP levels O/5/2, 1^3/2 ; 



lpi/2, and Ogg/2 located above the doubly-closed 8 Ca core. The aim of the present study is to 
compare the results of this calculation with those of Ref. [19], so as to try to clarify the role of 
the Ogg/2 orbital in heavy-mass calcium isotopes. 

The paper is organized as follows. In section 2 we give a brief outline of our calculations and 
some details about the choice of the SP energies. Section 3 is devoted to the presentation of the 
results for calcium isotopes beyond 48 Ca up to A=56. In section 4 we make some concluding 
remarks. 

2. Outline of calculations 

Our shell-model effective interaction has been derived within the framework of perturbation 
theory [17] starting from the CD-Bonn NN potential [22]. The high-momentum repulsive 
components of the bare NN potential have been renormalized by employing the so-called Vjow-k 
approach [20, 21], so as to obtain a smooth potential that preserves exactly the onshell properties 
of the original NN potential. Next, the TBME have been derived by way of the Q-box plus 
folded-diagram method [17], the Q-box being a collection of irreducible valence-linked Goldstone 
diagrams which we have calculated through third order in Vjow-k- 
The effective interaction V e s can be written in an operator form as 

V cS = Q-Q> J Q + Q> J Q J Q-Q> J Q J Q J Q+ ... , (1) 

where the integral sign represents a generalized folding operation, and Q' is obtained from Q by 
removing terms of first order in Vjow-k- The folded-diagram series is summed up to all orders 
using the Lee-Suzuki iteration method [23]. 

We have derived our shell-model effective interaction for the neutron-neutron, proton-proton, 
and proton- neutron channels using a model space spanned by the four neutron SP levels 0/ 5 / 2 , 
1^3/2) 1^1/2) Ogg/2, and by the sole proton SP level 0/ 7/ / 2 , all of them located above the doubly- 
closed 48 Ca core. The complete list of the matrix elements is reported in table 1. As already 
mentioned in the Introduction, in Ref. [19] the effective hamiltonian has been derived in a model 
space including the four pf orbitals outside doubly-closed 40 Ca, that is the standard choice of 
the model space in the study of calcium isotopes. 

As regards the choice of the neutron SP energies, we refer to [24] for those of the 0/ 5 / 2 , 
lP3/2> lPi/2 orbitals, obtained calculating the centroids of the observed single-neutron strengths 
in 49 Ca. Since there is no experimental evidence of a SP J n = | + state in 49 Ca, the 0<7g/2 SP 

energy has been fixed in order to reproduce the energy of the SP r = | state in 57 Ni [25]. 
It is worth mentioning that this choice leads also to a good agreement with experiment for the 
energies of J n = | + states in 51 Ti and 53 ~ 57 Cr. The adopted SP energies are e 3 / 2 = -4.60 MeV, 
e 1/2 = -2.86 MeV, e 5/2 = -1.20 MeV, and e 9/2 = -0.55 MeV. 

It should be pointed out that for protons the Coulomb force has been explicitly added to the 
Viow-k before deriving V c g. 

3. Results 

In this section, we report the results of the present calculations for the heavy-mass Ca isotopes 
from 50 Ca to 56 Ca (labelled as I) and compare them with those of [19] (labelled as II), where 
40 Ca was assumed as a closed core. The calculations have been performed by using the Oslo 
shell- model code [26]. 

The calculated and experimental [6] energies of 50 Ca are reported in table 2. We see that both 
calculations reproduce quite well the observed energies of positive-parity states. The negative- 
parity states can be obviously obtained only in calculation I, and the interpretation of their 
nature is strictly connected with the corresponding states in the 48 Ca energy spectrum. 



Table 1. Proton-proton, neutron- neutron, and proton- neutron matrix elements (in MeV). They 
are antisymmetrized and normalized. 



nbhjb n c l c j c n d ldjd J T z TBME 



0/7/2 0/7/2 0/7/2 0/7/2 





1 


-1.867 


0/7/2 0/7/2 0/7/2 0/7/2 


2 


1 


-0.483 


0/7/2 0/7/2 0/7/2 0/7/2 


4 


1 


0.160 


0/7/2 0/7/2 0/7/2 0/7/2 


6 


1 


0.479 


0/5/2 0/5/2 0/5/2 0/5/2 





-1 


-0.841 


O/5/2 O/5/2 lP3/2 l£»3/2 





-1 


-0.708 


O/5/2 O/5/2 lpi/2 lPl/2 





-1 


-0.396 


O/5/2 O/5/2 0^9/2 O39/2 





-1 


1.922 


1P3/2 1P3/2 1^3/2 1^3/2 





-1 


-0.802 


1^3/2 1P3/2 lPl/2 lPl/2 





-1 


-1.028 


!P3/2 !P3/2 0^9/2 0^9/2 





-1 


1.305 


lPl/2 lPl/2 lPl/2 lPl/2 





-1 


-0.161 


lPl/2 lPl/2 O59/2 O59/2 





-1 


0.980 


059/2 0#g/2 0^9/2 059/2 





-1 


-1.067 


O/5/2 1P3/2 O/5/2 1P3/2 


1 


-1 


0.138 


O/5/2 1^3/2 1P3/2 lPl/2 


1 


-1 


-0.061 


!P3/2 lPl/2 !P3/2 lPl/2 


1 


-1 


0.244 


0/5/2 O/5/2 O/5/2 O/s/2 


2 


-1 


-0.207 


O/5/2 O/5/2 O/5/2 lf> 3 /2 


2 


-1 


0.089 


O/5/2 O/5/2 O/5/2 lpi/2 


2 


-1 


-0.232 


O/5/2 O/5/2 1^3/2 1P3/2 


2 


-1 


-0.126 


O/5/2 O/5/2 lp 3 /2 lPl/2 


2 


-1 


-0.285 


O/5/2 O/5/2 0^9/2 O^g/2 


2 


-1 


0.266 


O/5/2 1P3/2 O/5/2 1P3/2 


2 


-1 


0.219 


O/5/2 1P3/2 O/5/2 lPl/2 


2 


-1 


0.353 


O/5/2 1P3/2 1P3/2 1P3/2 


2 


-1 


0.133 


0/5/2 1P3/2 lP3/2 lPl/2 


2 


-1 


0.221 


O/5/2 1P3/2 0^9/2 0^9/2 


2 


-1 


-0.554 


O/5/2 lPl/2 O/5/2 lPl/2 


2 


-1 


-0.197 


O/5/2 lPl/2 1P3/2 1P3/2 


2 


-1 


-0.127 


O/5/2 lPl/2 1P3/2 lPl/2 


2 


-1 


-0.386 


O/5/2 lPl/2 0^9/2 0^9/2 


2 


-1 


0.684 


!P3/2 !P3/2 !P3/2 !P3/2 


2 


-1 


-0.308 


J-P3/2 L P3/2 L PZ/2 L Pl/2 


9 


1 

_L 


\J.OO^± 


1P3/2 lP 3 /2 0^9/2 Ofifg/2 


2 


-1 


0.477 


!P3/2 lPl/2 1 P3/2 lPl/2 


2 


-1 


-0.452 


1P3/2 lPl/2 0^9/2 Ofifg/2 


2 


-1 


0.397 


Oc/9/2 0^9/2 059/2 0^9/2 


2 


-1 


-0.771 


O/5/2 1^3/2 O/5/2 1P3/2 


3 


-1 


0.163 


O/5/2 1P3/2 O/5/2 lPl/2 


3 


-1 


0.044 


O/5/2 lPl/2 O/5/2 lPl/2 


3 


-1 


0.307 


0/5/2 O/5/2 O/5/2 O/5/2 


4 


-1 


0.198 


O/5/2 O/5/2 O/5/2 lj> 3 /2 


4 


-1 


0.207 


O/5/2 O/5/2 0g 9/2 0g 9 /2 


4 


-1 


0.228 


0/5/2 1P3/2 O/5/2 1P3/2 


4 


-1 


-0.243 


O/5/2 1P3/2 0^9/2 Offg/2 


4 
4 


-1 
-1 


-0.613 
-fl 977 
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O/5/2 O59/2 O/5/2 0g 9/2 


2 


-1 


-0.569 


O/5/2 0^9/2 O/5/2 0^9/2 


3 


-1 


-0.436 


O/5/2 0^9/2 !P3/2 0^9/2 


3 


-1 


0.243 


1P3/2 0^9/2 1^3/2 Ofifg/2 


3 


-1 


-0.827 


O/5/2 0#g/ 2 O/5/2 0g 9/2 


4 


-1 


0.062 


O/5/2 O59/2 IP3/2 Ofifg/2 


4 


-1 


0.170 


O/5/2 O59/2 IP1/2 O^g/2 


4 


-1 


0.052 


1^3/2 Ofifg/2 1P3/2 Ofifg/2 


4 


-1 


0.060 


1P3/2 Ofifg/2 lPl/2 O^g/2 


4 


-1 


-0.059 


lPl/2 O^g/2 lPl/2 Ofifg/2 


4 


-1 


0.112 


0/5/2 0#g/2 O/5/2 Offg/2 


5 


-1 


-0.178 


O/5/2 0^9/2 1^3/2 Ofifg/2 


5 


-1 


0.196 


O/5/2 O59/2 lPi/2 O^g/2 


5 


-1 


-0.275 


!P3/2 O^g/2 !P3/2 Ofifg/2 


5 


-1 


-0.191 


1^3/2 05g/2 lPl/2 O^g/2 


5 


-1 


0.487 


lPl/2 Ofifg/2 lPl/2 O^g/2 


5 


-1 


-0.401 


O/5/2 0#g/2 O/5/2 0^9/2 


6 


-1 


0.153 


O/5/2 0^9/2 1^3/2 O59/2 


6 


-1 


0.194 


IP3/2 O^g/2 lp 3 /2 Ofifg/2 


6 


-1 


0.400 


0/5/2 O^g/2 O/5/2 Offg/2 


7 


-1 


-0.945 


O/7/2 O/5/2 O/7/2 O/5/2 


1 





-2.430 


O/7/2 O/5/2 O/7/2 O/5/2 


2 





-1.562 


O/7/2 O/5/2 O/7/2 lp 3/2 


2 





-0.638 


0/7/2 1P3/2 O/7/2 1P3/2 


2 





-0.830 


O/7/2 O/5/2 O/7/2 O/5/2 


3 





-0.422 


O/7/2 O/5/2 0/ 7 /2 lp 3/2 


3 





0.169 


O/7/2 O/5/2 O/7/2 lpi/2 


3 





-0.295 


O/7/2 1P3/2 O/7/2 lp 3 /2 


3 





-0.369 


O/7/2 1P3/2 0/7/2 lPl/2 


3 





0.790 


0/7/2 lPl/2 O/7/2 lPl/2 


3 





-0.678 


O/7/2 O/5/2 O/7/2 O/5/2 


4 





-0.891 


O/7/2 O/5/2 O/7/2 lp 3/2 


4 





-0.189 


O/7/2 O/5/2 O/7/2 lpi/2 


4 





-0.546 


u /7/2 J-P3/2 u /7/2 -LP3/2 


A 


n 

u 


u.uoo 


O/7/2 1P3/2 O/7/2 lPl/2 


4 





-0.344 


O/7/2 lPl/2 O/7/2 lPl/2 


4 





-0.573 


O/7/2 O/5/2 O/7/2 O/5/2 


5 





0.157 


O/7/2 O/5/2 O/7/2 lp 3 / 2 


5 





0.246 


0/7/2 1P3/2 O/7/2 lp 3 /2 


5 





-1.027 


O/7/2 O/5/2 O/7/2 O/5/2 


6 





-1.796 


O/7/2 Ofifg/2 O/7/2 0c/9/ 2 


1 





-1.138 


O/7/2 Ofifg/2 O/7/2 O59/2 


2 





-0.846 


O/7/2 O^g/2 O/7/2 0^9/2 


3 





-0.288 


O/7/2 Ofifg/2 O/7/2 Offg/2 


4 





-0.331 


O/7/2 O59/2 O/7/2 0g 9 /2 


5 





-0.066 


O/7/2 Ofifg/2 O/7/2 O59/2 


6 
7 






-0.499 
n nifi 



As a matter of fact, in 8 Ca the J n = 3f , 4^, 5j~ states, which can be interpreted as lp— Ih 
Z = 20 cross-shell excitations, are located ~ 0.5 MeV higher in energy than the corresponding 
ones in 50 Ca. This may be seen as an indication that, when increasing the number of neutrons, 
the structure of the above negative-parity states changes because of configuration mixing with 
states built on the neutron 0g 9 / 2 orbital. Our results seem to indicate that the overlap with 
our model space is larger for J 71 = 5j" than for J 71 = 3j" , 4^. The experimental J 77 = 2f, 4f 
states have no theoretical counterpart in calculation I since their structure is dominated by 
the ^(0/7/2) (IP3/2) 3 configuration [6, 19]. Note that the results for positive-parity states with 
calculations I and II are close to each other, thus confirming the soundness of the perturbation 
theory of the shell-model effective interaction. We have also calculated the B(E2;2f — > 0^) 
transition rate in 50 Ca employing effective operators obtained at third order in perturbation 
theory, consistently with the derivation of H e g as in Ref. [19]. The calculated value is 7.1 e 2 fm 4 
to be compared with the experimental one 7.5 ±0.2 e 2 fm 4 [27], while with calculation II we have 
obtained 10.9 e 2 fm 4 [19]. 

Table 2. Calculated energy levels (in MeV) of 50 Ca and 52 Ca compared with the corresponding 
experimental ones [6, 28] up to 5.2 MeV excitation energy. The energies of the states whose 
parity and/or angular momentum are uncertain are reported in parenthesis. 





Expt. 


50 Ca 
Calc. I 


Calc. II 




Expt. 


52 Ca 
Calc. I 


Calc. II 


0+ 


0.000 


0.000 


0.000 


0+ 


0.000 


0.000 


0.000 


2+ 


1.026 


0.939 


0.953 


2+ 


2.562 


2.358 


2.396 


2+ 


(3.004) 


2.856 


2.941 


1+ 


(3.150) 


3.205 


3.135 


1+ 


(3.532) 


3.472 


3.547 


3 - 


(3.992) 


4.702 




3 _ 


(3.998) 


4.698 












2+ 


(4.036) 




3.582 










0+ 


(4.470) 


4.839 


4.926 










4+ 


4.515 




4.567 










4 _ 


(4.832) 


5.590 












2+ 


(4.870) 


4.957 


5.231 










5" 


(5.110) 


5.188 













In table 3 we have reported the experimental and calculated energy spectra of Ca. We see 
that also for this nucleus the comparison between theory and experiment is remarkably good, 
apart from the J n = ^ state. The structure of this state is mainly obtained coupling the 
J w = 3i state in 50 Ca to a neutron in the 1^3/2 orbital, which explains the overestimation of 

the calculated energy by ~ 450 keV. In calculation I we find that the first excited J n = | + state 
is essentially a seniority v = 1 state. This lends further support to our choice of the 0g 9/ / 2 SP 

energy. It should be mentioned that in Ref. [8] both J 71 " = and J 71 " = are interpreted as 
lp — lh Z = 20 cross-shell excitations. 

As regards the negative-parity states, the agreement between calculations I and II is very 
good. Note that the J n = | state is outside the model space of calculation I since its structure 
is mainly made up by one neutron-hole in 48 Ca core and 4 neutrons filling the 1^3/2 orbital. 



The nucleus 52 Ca is the last one of this isotopic chain for which some experimental excited 
states have been observed [4-6]. The experimental energies of these states are reported in table 
2 together with the corresponding calculated ones up to 5.2 MeV. Both calculations I and II 
reproduce quite well the subshell closure of the lj>3/2 orbital. 

Table 3. Energy levels (in MeV) of 51 Ca and 53 Ca. The calculated 51 Ca energy levels are 
compared with the corresponding experimental ones [6,8]. The energies of the states whose 
parity and/or angular momentum are uncertain are reported in parenthesis. The 53 Ca energy 
levels are reported up to 4.0 MeV excitation energy. 







51 Ca 






53 Ca 






Expt. 


Calc. I 


Calc. II 




Calc. I 


Calc. II 


3~ 
2 


(0.000) 


0.000 


0.000 


1- 

2 


0.000 


0.000 


1- 
2 


(1.721) 


1.753 


1.586 


5~ 
2 


2.112 


2.039 


5~ 
2 


(2.379) 


2.190 


2.269 


3- 
1 + 


2.169 


2.402 


3~ 
2 


(2.937) 


2.829 


2.932 


2 


2.604 




7- 
2 


(3.437) 




3.429 








5~ 
2 


(3.479) 


3.564 


3.552 








7 + 


(3.845) 


4.309 










i+ 

2 


(4.155) 


4.126 










9~ 
2 


(4.322) 


4.281 


4.435 









For the heavier Ca isotopes, 53 ' 54 ' 55 Ca, only the ground states have been identified with spin 
and parity assignments [7]. In tables 3 and 4 we have reported our predicted spectra for these 
three nuclei and for 56 Ca up to 4.0, 4.5, 2.4, and 3.0 MeV excitation energy, respectively. From 
the inspection of the above tables, it is clear that the disagreement between calculations I and 
II increases with mass number, since the role of orbital becomes more and more relevant 
when adding valence neutrons. This has also repercussions on the two-neutron separation energy 
S2n for 56 Ca which is about 2 MeV larger in calculation I, as can be inferred from the inspection 
of figure 1 where the calculated ground-state energies of the even-mass isotopes relative to 48 Ca 
are compared with the measured (A = 50, 52) and estimated ones (A = 54, 56) [29]. 



Table 4. Theoretical energy levels (in MeV) of 54 Ca, 55 Ca, and 56 Ca. See text for details. 



r 


54 Ca 
Calc. I 


Calc. II 




55 Ca 
Calc. I 


Calc. II 


r 


56 Ca 
Calc. I 


Calc. II 


0+ 


0.000 


0.000 


5- 


0.000 


0.000 


0+ 


0.000 


0.000 


2+ 


2.587 


2.061 


2 


0.208 




2+ 


2.268 


0.963 


0+ 


2.783 


3.138 


\- 
2 


0.237 


1.104 


7- 


2.514 




5- 


3.116 




3~ 


1.924 


1.541 


2- 


2.745 




3+ 


3.157 


2.748 


1 + 


2.131 




3- 


2.825 




4- 


3.430 




5+ 
2 


2.204 




5- 


2.960 




2+ 


4.517 


4.172 


7- 
2 


2.325 


1.931 


4+ 


2.961 


1.487 








7+ 
2 


2.368 




2+ 


2.982 


2.224 








5- 
2 


2.381 


2.390 










4. Concluding remarks 

In summary, we have given here a shell-model description of heavy calcium isotopes assuming 
Ca as a closed core and including the neutron 0gg/ 2 orbital. This has been done by 
renormalizing the NN potential CD-Bonn by way of the Vj ow _k approach and then deriving 
the shell-model hamiltonian within the framework of perturbation theory. We have compared 
the results of the present study with those of the previous one [19] which started from doubly- 
closed 40 Ca and employed the lpOf model space. 

The results of both calculations are in good agreement with the available experimental data, 
and do not differ significantly up to 53 Ca. Actually, from 54 Ca on the differences between the 
energy spectra calculated within the two different model spaces increase with the number of 
valence neutrons. Clearly, new experimental data are needed to establish the real role of the 
neutron O59/2 orbital in the description of the structure of nuclei north-east of 48 Ca core. 
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